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We review recent results on QCD at finite temperature and non-vanishing baryon number density. We focus 
on observables which are of immediate interest to experimental searches for the Quark Gluon Plasma, i.e. the 
phase transition temperature, the equation of state in two and three flavour QCD and thermal effects on hadron 
masses. Some new developments in finite density QCD are also presented. 



1. Introduction 

1.1. Preface 

There are many reasons to study a complicated 
interacting field theory like QCD under extreme 
conditions, eg. at finite temperature and/or non- 
vanishing baryon number density. Of course, first 
of all we learn about the collective behaviour 
of strongly interacting matter, its critical be- 
haviour, equation of state and thermal properties 
of hadrons as well as their fate at high tempera- 
ture. However, in doing so we also have a look at 
the complicated non-perturbative structure of the 
QCD vacuum from a different perspective and can 
learn about the mechanisms behind confinement 
and chiral symmetry breaking. This does allow 
to check concepts that have been developed to 
explain these non-perturbative features of QCD. 

This review focuses on the former aspect of fi- 
nite temperature QCD studies; the latter has re- 
cently been discussed in In fact, an impor- 
tant step in the direction of understanding con- 
finement in terms of the dual superconductor pic- 
ture has been taken recently by the Pisa group . 
They analysed the scaling behaviour of an order 
parameter for monopole condensation and could 
show that it scales in SU(2) and SU(3) gauge the- 
ories with the critical exponents expected for the 
deconfinement transition (Fig. This is a strong 
hint for the survival of the monopole condensa- 
tion mechanism as a source for confinement in 
the continuum limit. 



*The work has been supported by the TMR network 
ERBFMRX-CT-970122 and the DFG under grant Ka 
1198/4-1 . 
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Figure 1. Finite size scaling analysis of an order 
parameter for monopole condensation at the de- 
confinement transition temperature of the SU(2) 
gauge theory ||,|. 

1.2. Outline 

Significant progress in finite temperature lat- 
tice QCD is based on the use of improved ac- 
tions. In the pure gauge sector we have seen 
that these actions strongly reduce finite cut-off ef- 
fects and allow to calculate thermal quantities on 
rather coarse lattices in good agreement with con- 
tinuum extrapolated results. A similar approach 
is now followed in calculations with light quarks. 
The use of improved staggered and Wilson actions 
leads to new estimates of the critical temperature 
and the equation of state which we will discuss in 
the following sections. Moreover, thermodynamic 
calculations with domain wall fermions (DWF), 
which can greatly improve the chiral properties 
of lattice fermion actions, reached a stage where 
first quantitative results can be reported on 
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and thermal masses. The latter will be discussed 
in Section 4. In Section 5 we will discuss new de- 
velopments in the analysis of QCD at finite den- 
sity. Section 6 contains an outlook. 

An important part of finite temperature (T) 
and density {ub) (chemical potential (/i)) calcu- 
lations consists of trying to understand the QCD 
phase diagram, i.e. the order of the phase tran- 
sition as function of the number of flavours (ny) 
and quark masses (m^) and its dependence on T 
and Ub (or fi). We will not discuss these issues 
here. The finite temperature phase diagram has 
been discussed extensively in recent reviews |^ 
and the possibly rather complex phase structure 
of QCD at high density and low temperatures has 
been the topic of E. Shuryak contribution to this 
conference 

2. The Critical Temperature 

One of the basic goals of lattice QCD calcula- 
tions at finite T is to provide quantitative results 
for the phase transition temperature. In the pure 
gauge sector this goal has been achieved. The 
value of the critical temperature for the decon- 
fining phase transition in a SU{3) gauge theory 
is known with small statistical errors. Remaining 
systematic errors are of the order of 3%. They 
arise from the calculation of the string tension at 
T = which is used to set the scale for Tc. In 
Table 1 we summarise results for Tc/ ^/a obtained 
with different gauge actions on lattices with vary- 
ing temporal extent ~ 1/ (aT) and extrapo- 
lated to the continuum limit (iV^ oo). 

Using for the string tension the value^ -^/ct ~ 

^This may be deduced from quenched spectrum calcula- 

Table 1 

Critical temperature of the SU{3) gauge theory 
in units of ^/a obtained from calculations with 
different gauge actions and extrapolated to the 
continuum limit. 



action 


Ref. 




standard Wilson 
Symanzik impr. (tree level) 
RG-improved 


[5,6] 

[6] 

[7] 


0.630 (5) 
0.634 (4) 
0.650 (5) 



425 MeV, we find ~ 270 MeV. This large 
value can be understood in terms of the particle 
spectrum of the quenched theory; in the low tem- 
perature, confining phase there exist only rather 
heavy glueballs (tog>1.7 GeV). A rather large 
temperature thus is needed to create a sufficiently 
dense glueball gas, which can trigger a deconfin- 
ing transition. 

Investigations of the critical temperature in 
QCD with quarks of finite mass indeed have 
shown that the transition temperature drops 
rapidly with decreasing quark masses]^ 

Unlike in the pure gauge theory the transition 
temperature for QCD with finite quark masses 
does seem to be strongly affected by the dis- 
cretization scheme used for the fermion action. 
The early calculations with the standard stag- 
gered Q and Wilson [|I^ actions led to widely 
different estimates for Tc- Finite-cut off effects 
are thus expected to be large and improvement 
of the fermion action should be expected to be 
important. Indeed, the new calculations per- 
formed with improved Wilson (Clover) [p^ , p"4|jr^ ] 
and staggered [|l^ actions as well as with DWF 
[|l7|| yield transition temperatures which are in 
much better agreement among each other. Unfor- 
tunately, this statement is, at present, only par- 
tially correct. In fact, for comparable values of 
TTiq, i.e. fixed ratios of pseudo-scalar ("pion") and 
vector meson ("rho") masses, mps/mv, it only 
holds when we set the scale for Tc using a hadron 
mass, eg. my- When we follow a similar ap- 
proach and use the string tension, y/a, to set the 
scale the agreement is less evident. In this case 
calculations with the Wilson fermion action typ- 
ically lead to transition temperatures about 20% 
below the results obtained with staggered fermion 



tions {rnpl\f(j = 1.81 (4) |8|). Recent estimates from CP- 
PACS data seem to lead to a somewhat smaller value |b|. 
^At intermediate values of rriq the transition only reflects 
a rapid cross-over in thermodynamic observables rather 
than a true phase transition. The peak in the chiral sus- 
ceptibility or the Polyakov-Loop susceptibility defines in 
this case a pseudo-critical temperature. We will continue 
to call this temperature the transition temperature at fi- 
nite values of niq For 2-fiavour QCD it is expected to be a 
critical temperature, corresponding to a phase transition, 
only in the chiral limit. For 3-flavour QCD it is a c r itical 
temperature below a certain critical quark mass |U)lhl|. 
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Figure 2. Transition temperatures for 2-flavour QCD in units of mv (left) and ^/a (right), respectively. 
The corresponding quark masses for staggered fermions are indicated in the right frame. Filled (open) 
symbols are for results from Nr — 4 (6) lattices. The large dots drawn for mps /mv = are just indicative 
for the chiral limit. They would correspond to a critical temperature of Tc — 170 MeV. The long vertical 
lines indicate the location of the physical limit, raps = rriTr — 140 MeV. The solid lines connect a data 
set obtained with the standard staggered ferniion action ||l8|. They are meant to guide the eye. 



actions. The current status of the determination 
of Tc for 2-flavour QCD is summarised in Fig. |2[ 
As can be seen the simulations with Clover 
fermions performed by different groups are con- 
sistent with each other. Results on Tc/my do 
not seem to depend significantly on the gauge 
action (one plaquette Wilson p4|, Symanzik im- 
proved jl^] or RG-improved |15[ ) nor does it 
seem to be important whether tadpole p^ , p^ 
or non-perturbative Clover coefficients are 
used. Also the DWF results |jl^ seem to be in- 
sensitive to the gauge action chosen (plaquette or 
RG-improved) and are consistent with results ob- 
tained with the Clover action. Results obtained 
with an improved staggered fermion action [p^ , 
the p4-actionQ also agree with these data within 
10%. 

Unfortunately this consistent picture is, at 
present, not reproduced when calculating Tc in 
units of y/a (right frame of Fig. ||). A possi- 
ble source for this discrepancy might be the cal- 
culation of the heavy quark potential, which in 
the case of Wilson fermions so far has only been 
performed on rather small spatial lattices, e.g. 

*The p4-action improves the rotational symmetry of the 
quark propagator in C'(p'*). It strongly reduces the cut-off 
effects in bulk thermodynamic observables ted], see Fig. 4. 



8"^ X 16. This may lead to an overestimate of 
the string tension. It, however, also is possible 
that calculations on Nr = 4 lattices are still per- 
formed at too strong coupling and do not allow 
for a unique determination of the scale. A hint 
in this direction may be the large difference in 
Tcl\f^ observed in calculations on TV,- = 4 and 
6 lattices with the non-perturbatively improved 
Clover action [^. Clearly more work is needed 
here to establish a unique result for Tc using dif- 
ferent fermion formulations. 

In general, we note that the transition temper- 
ature obtained with improved actions tends to be 
larger than what previously has been quoted on 
the basis of calculations performed with the stan- 
dard staggered action. If the quark mass depen- 
dence does not change drastically closer to the 
chiral limitj^ the current data suggest 

Tc ~ (170- 190) MeV (1) 

for 2-flavour QCD in the chiral limit. In fact, 
this estimate also holds for 3-flavour QCD. Calcu- 
lations which are currently performed for Uf = 2 
and 3 using the same improved staggered fermion 

^Calculations within the framework of quark- rucaon mod- 
els suggest a rapid drop of Tc for rnpg < |19| . 
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Figure 3. The transition temperature in 2 
(square) and 3 (circles) flavour QCD versus 
'mpsl\fo using the p4-action. The dashed 
band indicates the uncertainty on Td \fo in the 
quenched limit. 

action, suggest that the flavour dependence of 
is rather weak (see Fig. ||). 

It is remarkable that the transition tempera- 
ture drops significantly already in a region where 
all hadron masses are quite large. This is appar- 
ent from Fig. ^ where we show in units of \fG 
plotted vs. mpsi \fo for Uf — 2 and 3. As can be 
seen the transition temperature starts deviating 
from the quenched values for ?nf s<(6 — 7)^/ct ~ 
2.5 GeV. We also note that the dependence of Tc 
onmps/ \fa is almost linear in the entire mass in- 
terval. This might be expected for light quarks in 
the vicinity of a 2"'' order chiral transition where 
the pseudo critical temperature depends on the 
mass of the Goldstone-particle like 



Tcim^) - Tc(0) ~ TO 



(2) 



For 2-flavour QCD where the critical indices are 
expected to belong to the universality class of 3- 
d, 0(4) symmetric spin models one would indeed 
expect 1/(36 = 1.1. However, this clearly cannot 
be the origin of the quasi linear behaviour ob- 
served for large hadron masses independent^ of 
Hf. A resonance gas model would probably be 
more appropriate to describe the thermodynam- 



ics for these heavy quarks. 

3. The Equation of State 

In the pure gauge sector bulk thermodynamic 
quantities such as the pressure and energy density 
have been analysed in detail. It has been verified 
that the most significant cut-off effects result from 
high momentum modes, which dominate the in- 
finite temperature, ideal gas limit. Improved ac- 
tions that lead to small cut-off effects in the ideal 
gas limit]] still do so at finite temperature 

In a recent analysis the CP-PACS collaboration 
calculated the pressure and energy density using 
the RG-improved (Iwasaki) action They con- 
firm that after extrapolation to the continuum 
limit also this action yields results for the SU (3) 
equation of state which are, within an error of 
(3-4)%, consistent with the continuum extrapo- 
lation obtained from the Wilson |2^], improved 
as well as some fixed point actions |2^. In 
fact, the Wilson and RG-improved actions repre- 
sent extreme cases for such calculations. Their 
cut-off corrections have opposite sign which on 
coarse {Nr = 4) lattices leads to deviations of the 
pressure by about 25% from the continuum ex- 
trapolated result. In view of this the agreement 
reached with different discretization schemes is 
rather reassuring. 

Also in the presence of light quarks the use 
of an improved action thus seems to be manda- 
tory, if one wants to calculate bulk thermody- 
namic quantities. The standard staggered and 
Wilson actions are known to lead to large devi- 
ations from the continuum ideal gas behaviour 
on coarse lattices (small N^) For staggered 
fermions improved actions, leading to smaller cut- 
off effects in thermodynamic observables, can be 
constructed by adding suitably chosen three-link 
terms to the conventional one- link terms [^2| , p5| . 
In Fig. ^ (left frame) we show the cut-off de- 
pendence of the ideal gas pressure obtained from 
these minimally improved staggered fermion ac- 
tions (Naik mi and p4 pi action)|. 



"A similar conclusion holds for nf = 2, when one analyses 
the unimproved standard staggered fermion data. 



'^In the ideal gas lij 

ysed analytically [ 
fermion actions see 

* Staggered actions with even smaller cut-off effects have 



the cut-off dependence can be anal- 
For a recent analysis of staggered 
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Figure 4. Cut-ofF dependence of the pressure of an ideal Fermi gas calculated on spatially infinite lattices 
with temporal extent and different staggered fermion actions (left) and the pressure for 2-flavour QCD 
calculated by using the standard staggered fermion action ||2^] and the p4-improved action (right). 
The dashed band shows an estimate for the 2-flavour pressure in the continuum limit (see text). The 
horizontal line for T < Tc shows the pressure of a massless pion gas. 



As expected the systematics seen for the cut-off 
dependence in the ideal gas limit carries over to fi- 
nite T. In the right frame of Fig. ^ we show results 
for the pressure of 2-flavour QCD obtained with 
the standard staggered action (and Wilson gauge 
action) on lattices with temporal extent Nr = 4 
and 6 These are compared with results ob- 
tained with the p4-action (and a Symanzik im- 
proved gauge action) [|l]j2^. In the latter case 
also fat links ||2^ have been introduced in the 
one-link terms to improve the flavour symmetry 
of the staggered action. 

While the pressure calculated with the stan- 
dard staggered action rapidly overshoots the ideal 
gas limit (as expected from the analysis of the 
cut-off effects in the ideal gas limit) the results 
obtained with the p4-action stay below the ideal 
gas limit and show a temperature dependence 
very similar to what has been found in the pure 
gauge sector. We note that the calculations with 
the p4-action have been performed with rather 
large quark masses, m/T = 0.4, corresponding to 
nips/mv — 0.7, while the calculation with stan- 
dard staggered action are for mps/my — 0.3. In 
the high temperature phase this does not seem 



been constructed |26| , p7| ] . However, they also require a 
significantly larger numerical effort. 



to constitute a major problemQ Smaller quark 
masses are, however, definitely needed close to 
Tc and below in order to become sensitive to the 
contributions from light pion modes. 

In the case of the SU(3) gauge theory the mag- 
nitude of the cut-off dependence in the tempera- 
ture range up to a few times Tc has been found to 
be about half of what has been calculated for the 
ideal gas, infinite temperature limit. If this con- 
tinues to hold for the fermion sector one should 
expect, that the current result, obtained with 
the p4-action on = 4 lattices, underestimates 
the continuum result by about 15% for T>Tc. 
Based on this consideration an estimate for the 
continuum extrapolated pressure is also shown in 
Fig. |4| We also note that a calculation performed 
with identical, improved staggered fermion ac- 
tions shows that the pressure normalised to the 
continuum ideal gas value has the same temper- 
ature dependence in 2 and 3 flavour QCD p^ . 

Unfortunately, Wilson actions with similarly 
good high temperature behaviour have not been 
constructed so far. The Clover action does not 



^Calculations for 4-flavour QCD with quark masses 
m/T = n.2 and 0.4 show no significant quark mass depen- 
dence [B0||. Moreover, also in the continuum the pressure 
of an ideal gas of fermions with mass m/T = 0.4 deviates 
by less than 10% from that of a massless gas. 
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Figure 5. Pressure vs. T/Tc calculated with Wil- 
son fermions on a IQ^ x 4 lattice for two different 
meson mass ratios, raps/mv- 

improve the ideal gas behaviour, i.e. it has the 
same infinite temperature limit as the Wilson ac- 
tion. Nonetheless, a first attempt to calculate 
the equation of state, using the tadpole-improved 
Clover action combined with the RG-improved 
gauge action, has now been undertaken by the 
CP-PACS collaboration A preliminary re- 

sult for the pressure in 2-flavour QCD is shown in 
Fig. 1^. Similar to simulations with the standard 
staggered fermion action one observes an over- 
shooting of the ideal gas limit reflecting the cut- 
off effects in the unimproved fermion sector. 

4. Thermal Masses 

Understanding the temperature dependence of 
hadron properties, e.g. their masses and widths, 
is of central importance for the interpretation of 
heavy ion experiments. Thermal modifications of 
the heavy quark potential influence the spectrum 
of heavy quark bound states. Their experimen- 
tally observed suppression |^l[ thus is expected to 
be closely linked to the deconfining properties of 
QCD above Tc |32|. Changes in the chiral conden- 
sate, on the other hand, influence the light hadron 
spectrum and may leave experimental signatures, 
for instance in the enhanced dilepton production 
observed in heavy ion experiments [p3[. 



In numerical calculations on Euclidean lat- 
tices one has access to thermal Green's functions 
G'/f(r, r) in fixed quantum number channels, 77, 
to which in particular at high temperature many 
excited states contribute. As the temporal di- 
rection of the Euclidean lattice is rather short at 
finite temperature one usually has not enough in- 
formation on the correlation function to extract 
reliably thermal effects on the (pole) masses. A 
way out may come from the use of anisotropic 
lattices, which has extensively been explored by 
the QCDTARO group ||||. In their recent anal- 
ysis with quenched Wilson fermions they show 
evidence for a rapid change in the thermal corre- 
lation functions across Tc which indicates a rapid 
change in thermal masses above T,.- At the same 
time their analysis of pseudo-scalar wavefunctions 
does, however, show evidence that the correlation 
among quarks in this quantum number channel 
remains strong. In how far this hints at the pres- 
ence of light " pionic" bound states has to be anal- 
ysed further. 

Anisotropic lattices in combination with 
NRQCD have also been used to analyse thermal 
effects on heavy quark bound states [|5|. Large 
mass shifts have been observed for the first ex- 
cited states. 

The missing information on the long distance 
behaviour of thermal correlation functions may 
also be overcome by using refined techniques to 
analyse the numerical data on thermal correla- 
tion functions. It recently has been suggested 
that the maximum entropy method may help 
to extract more detailed information on thermal 
modifications of hadronic spectral functions |Q . 

So far the most convincing evidence for modifi- 
cations of hadron properties at high temperature 
comes from the analysis of the behaviour of corre- 
lation function at large spatial separations, which 
yields screening masses [Q, as well as thermal 
susceptibilities, 

XH = dT j d\GH{r,r) , (3) 

with Gh denoting the hadron correlation func- 
tion in the quantum number channel H . If there 
is only a single stable particle of mass mn con- 
tributing to Gh then xh ^ ''^h^- general 
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Figure 6. Susceptibilities in different quantum 
number channels for 2-flavour QCD with stag- 
gered fermions on lattices of size 8^ x 4. The 
various points are obtained from a Ferrenberg- 
Swendsen reweighting analysis based on measure- 
ments at a few /3-values in the interval shown. 
The horizontal line indicates the location of /3c- 

the susceptibilities, however, define only effective 
masses, i.e. they average over contributions from 
the ground state and all excited states in a par- 
ticular quantum number channel. 

In Fig. H we show 1 /y^XH for 2-flavour QCD ob- 
tained with staggered fermions of mass = 0.02 
on lattices of size 8'^ x 4.. This figure is based 
on data from |39|. As can be seen the fo and 
TT masses become (almost) degenerate at Pc while 
the ao remains heavy at Pc- The former behaviour 
is expected, /o and tt correlation functions are re- 
lated through a rotation in SU{2) flavour space. 
The degeneracy thus reflects the restoration of 
the SU{2) flavour symmetry. The difference of 
Xao ^-iid Xtt on the other hand reflects the persis- 
tence of the [^4(1) symmetry breaking. A crucial 
question here is how the gap observed for non- 
zero TUq changes in the chiral limit. At high tem- 
perature topologically non-trivial gauge field con- 
figurations are expected to be suppressed, which 
in turn would lead to a strong reduction in the 
strength of [/^(l) symmetry breaking and thus 
in a strong reduction of the mass splitting be- 
tween ao and tt. For 2-flavour QCD it is expected 
that the quark mass dependence is quadratic. 



{mao-m^) ix-^-Xao) ~ A+Bml. Previous in- 
vestigations of this were, however, not conclusive 
1^. If a quadratic ansatz is assumed calculations 
with staggered fermions led to the conclusion that 
a non-zero mass splitting remains also above Tc, 
i.e. the f/A(l) remains broken ||]. The problem 
has now been addressed again by the Columbia 
group [0 using DWF. Due to the improved chi- 
ral properties of this action one should find a 
quadratic dependence on the quark mass. For 
T ~ 1.2Tc the Columbia group indeed does ob- 
serve such a quark mass dependence. In the chi- 
ral limit they find a non-zero mass splitting from 
the susceptibilities as well as from the analysis of 
screening masses [T^ ], 

- TO^ = 0.0606 (67) + 9.66 (58)to^ . (4) 

The J/a(1) thus remains broken above Tc, al- 
though the mass splitting is strongly reduced (see 
Fig. |6|). This picture is also supported by an 
analysis of the disconnected parts of flavour sin- 
glet correlation functions in quenched QCD using 
DWF . A vanishing of these would signal a de- 
generacy between the a and 77 mesons. The dis- 
connected parts get non-zero contributions only 
from topologically non-trivial configurations . 
Indeed, these are strongly suppressed above Tc. 
However, at T = 1.25rc about 10% of the con- 
figurations do still carry a non-zero topological 
charge indicating the persistence of Ua{^) sym- 
metry breaking above Tc. 

5. Finite Density QCD 

Finite density calculations in QCD are affected 
by the well known sign problem, i.e. the fermion 
determinant becomes complex for non-zero values 
of the chemical potential fi and thus prohibits the 
use of conventional numerical algorithms. The 
most detailed studies have been performed so far 
using the Glasgow algorithm , which is based 
on a fugacity expansion of the grand canonical 
partition function at non-zero fi. 



-'GC 



ifi/T,T, V) = 



aV 

E 

B=-aV 



z'^ZbIT, V) 



(5) 



where z — exp {fJ-/T) is the fugacity and Zb are 
the canonical partition functions for fixed quark 
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Figure 7. Polyakov loop expectation value (left) calculated on x 2 and the heavy quark potential 
(right) calculated on 16"^ x 4 lattices in quenched QCD at zero and non-zero baryon number, B/3. 



number B; a = 3,6 for one species of staggered 
or Wilson fermions, respectively. However, this 
approach so far did not overcome the severe nu- 
merical difficulties. 

It thus seems necessary to approach the finite 
density problems from another point of view. A 
reformulation of the original ansatz may lead to a 
representation of the partition function which, in 
the ideal case, would require the averaging over 
configurations with strictly positive weights only, 
or at least would lead to a strong reduction of 
configurations with negative weights. 

A big step away from the original formulation is 
to start from a Hamiltonian approach . Here 
it has been shown that problems with a fluctuat- 
ing integrand can successfully be reformulated in 
terms of a model where the configurations gener- 
ated do have strictly positive weights (meron clus- 
ter algorithm Q]). Whether such an approach 
can be applied to QCD remains to be seen. 

An alternative formulation of finite density 
QCD is given in terms of canonical rather 
than grand canonical partition functions ]45| , 
i.e. rather than introducing a non-zero chemi- 
cal potential through which the number density 
is controlled one introduces directly a non-zero 
baryon number (or quark number B) from which 
the baryon number density on lattices of size 
X Nr is obtained as ub/T^ = f{Nr/N^f. 

After introducing a complex chemical potential 
in Zgc (Eq. ||) the canonical partition functions 



can be obtained via a Fourier transform^, 
Zb{T, V) = ^ r d0 e^*^/^ Zgc(*0, T, V) .(6) 

Also this formulation is by no means easy to use 
in general, i.e. for QCD with light quarks. In 
particular, it also still suffers from a sign problem. 
It, however, leads to a quite natural and useful 
formulation of the quenched limit of QCD at non- 
zero density p6| . 

5.1. Quenched limit of finite density QCD 

It had been noticed early that the straightfor- 
ward replacement of the fermion determinant by 
a constant does not lead to a meaningful static 
limit of QCD |^^. In fact, this simple replace- 
ment corresponds to the static limit of another 
theory with an equal number of fermion flavours 
carrying baryon number B and —B, respectively 
poll . This should not be too surprising. When one 
starts with and takes the limit of inflnitely heavy 
quarks something should be left over from the de- 
terminant that represents the objects that carry 
the baryon number. In the canonical formulation 
this becomes obvious. For — )■ oo one ends 
up with a partition function, which for baryon 
number B /3 still includes the sum over products 
of B Polyakov loops, i.e. the static quark propa- 
gators which carry the baryon number . This 

^"The use of this ansatz for the calculation of canonical 
partition function s as cx pa,nsion coefficients for Zqc has 
been discussed in llTlkS . 
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limit also has some analogy in the grand canonical 
formulation where the coupled limit ruq , — s- oo 
with exp (/i)/2mq kept fixed has been performed 

miB. 

In the confined phase of QCD the baryon num- 
ber is carried by the rather heavy nucleons. Ap- 
proximating them by static objects may thus be 
quite reasonable and we may expect to get valu- 
able insight into the thermodynamics of QCD 
at non-zero baryon number density already from 
quenched QCD. 

From a numerical point of view there is hope 
that we can get control over this limit using dif- 
ferent approaches. At non-zero ^ a variant of the 
Glasgow approach Q seems to become applica- 
ble for large quark masses and also the static limit 
of Blum et al. may be explored further |Q . 

In the canonical approach simulations at non- 
zero B can be performed on relatively large lat- 
tices and the use of baryon number densities up 
to a few times nuclear matter density may be 
possible H^JI^. The simulations performed so 
far show the basic features expected at non-zero 
density. As can be seen from the behaviour of the 
Polyakov loop expectation value shown in Fig. |^ 
the transition region gets shifted to smaller tem- 
peratures (smaller coupling [3). The broadening 
of the transition region may suggest a smooth 
crossover behaviour at non-zero density. How- 
ever, in a canonical simulation it also may indi- 
cate the presence of a region of coexisting phases 
and thus would signal the existence of a 1** order 
phase transition. This deserves further analysis. 

Even more interesting is the behaviour of the 
heavy quark potential in the low temperature 
phase. As shown in the right frame of Fig. |^ 
the potential does get screened at non-vanishing 
number density. This will have a direct influence 
on heavy quark bound states at high density. 



'^^This is a well known limit in statistical physics. When 
deriving the non-relativistic gas limit from a relativistic 
gas of particles with mass m, the rest mass is splitted off 
from the chemical potential, = /inr + "i, in order to 
cancel the corresponding rest mass term in the particle 
energies. On the lattice m = ln(2m5) for large bare quark 
masses. 



6. Outlook 

We have focused in this review on the calcula- 
tion of basic thermodynamic quantities which are 
of immediate interest to experimental searches for 
the Quark Gluon Plasma. 

In quenched QCD the critical temperature and 
the equation of state have been calculated on the 
lattice and extrapolated to the continuum with 
an accuracy of a few percent. These calculations 
set a benchmark for many analytical studies of 
QCD thermodynamics ||5^. The progress made 
in developing and testing improved fermion ac- 
tions for thermodynamic calculations shows that 
a similar accuracy for QCD with light quarks is 
within reach. The current systematic studies with 
different improved fermion actions may soon lead 
to a determination of the transition temperature 
and the equation of state with similar accuracy. 

We do have reached some understanding of 
thermal effects on hadron properties. In partic- 
ular, modifications of the light meson spectrum 
due to flavour and approximate Ua{^) symme- 
try restoration have been established. However, 
lattice calculations so far did not come up with 
detailed quantitative results on thermal masses, 
which could be confronted with experimental 
data. There are a few promising ansatze which 
can lead to more detailed information on thermal 
modifications of hadronic spectral densities. 

Of course, there are many more important is- 
sues which have to be addressed in the future. 
Even at vanishing baryon number density we do 
not yet have a satisfactory understanding of the 
critical behaviour of 2-flavour QCD in the chi- 
ral limit and the physically realized situation of 
QCD with two light, nearly massless quarks and 
a heavier strange quark has barely been analysed. 

Moreover, the entire phase diagram at non-zero 
baryon number density is largely unexplored. An 
interesting phase structure is predicted at high 
density and low temperatures which currently is 
not accessible to lattice calculations. This does 
require new algorithmic developments. 

There are thus many interesting questions wait- 
ing to be answered in the next millennium. 
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